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STOKES' THEOREM FOR NONSMOOTH CHAINS 
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Abstract. Much of the vast literature on the integral during the last two centuries 
concerns extending the class of integrable functions. In contrast, our viewpoint is 
akin to that taken by Hassler Whitney [Geometric integration theory, Princeton Univ. 
Press, Princeton, NJ, 1957] and by geometric measure theorists because we extend the 
class of integrable domains. Let uj be an n-form defined on R m . We show that if cu is 
sufficiently smooth, it may be integrated over sufficiently controlled, but nonsmooth, 
domains 7. The smoother is oj, the rougher may be 7. Allowable domains include 
a large class of nonsmooth chains and topological n- manifolds immersed in R m . 
We show that our integral extends the Lebesgue integral and satisfies a generalized 
Stokes' theorem. 



requires both a smooth n-manifold M and a smooth (n— l)-form u). It was realized 
at some point that one side of Stokes' formula could be used to define the other side 
in more general situations. In particular Whitney [14] used the right side to define 
the left in some cases where the interior of M is smooth, even though the boundary 
DM is not smooth. More generally the interior could be piecewise smooth or a 
suitable limit of simplicial chains. Whitney systematically developed this insight 
by defining his flat norm on chains, based on rectilinear subdivisions of simplices 
of one higher dimension of which the chain in question is a partial boundary. He 
treated forms as cochains — linear functions on the vector space of chains. 

Stokes' theorem was thus extended by Whitney to integration of smooth forms 
over objects that were limits in the flat norm of chains. These include certain kinds 
of fractals. But other fractals escape Whitney's construction. Whitney's example 
of a function nonconstant on a connected arc of critical points [15] shows the limits 
to any generalization of Stokes' theorem. Stokes' theorem (which for arcs is just 
the fundamental theorem of calculus) must fail for such arcs. There exists a Jordan 
curve in 3-space that has Hausdorff dimension > 2, is not contained in any surface 
of finite area, and is not a limit in the flat norm of simplicial 1-chains. Whitney's 
methods to not define integration of forms over such a curve, nor do Lebesgue's. 
The methods for this paper accomplish exactly this. 
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I. Introduction 



The standard version of Stokes' theorem: 
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In a recursive construction we define a family of norms: For each real A with 
n < A < m the \-natural norm |^4|^ is defined on simplicial n-chains A in m-space. 
The definition utilizes projections of chains on hyperplanes, Whitney's method 
of subdivisions and partial boundaries, and Ath powers of the n-masses of the 
projected chains. The recursion is on the integer part of A — n. 

The following theorem implies that this is indeed a norm on oriented simplicial 
chains in M m (which is not obvious). Completion yields Banach spaces X nt \ = 
A"„ : A(K m )- These spaces are nested, containing increasingly complex chains as A 
increases. The parameter A acts as a dimension; it bounds the fractal complexity 
of elements of X 7lt \. The boundary operator d on simplices extends to a continuous 
operator d: X n ,\ — > X n - lyX - 

Theorem 1. Integration of an n-form u) defines a bounded linear functional L u 
on X 7lt \ if uj is of class C x ~ n . Furthermore, there is a constant c > such that for 
a simplicial n-chain A: 



[ lu < c |a|5jmu- 

Ja 



where \\u)\\\— n denotes the A — n norm on forms. 
This theorem extends results in [14] and [7]. 

Denote the dual space of cochains by X n ' X and the product of a cochain X E 
X n ' x with a chain A E X, hX by X • A. Define fT< A = {n-forms to of class C x ~ n 
defined on W n }. According to Theorem 1 the linear operator Ah J ui on simplicial 
n-chains in R m is bounded and determines an element L w E X n ' X . This results in 
a linear mapping L: fi" ,A — ► X n ' X . For A E X n ,\ wc define j A ui — • A. Special 
cases: If A is a simplicial n-chain (finite) in M m , then f A ui is identical to the 
Lebesgue integral for simplicial n-chains. If A is an n-chain (possibly infinite) in 
R", then J A u is identical to the Lebesgue integral for infinite n-chains with finite 
n-mass. Since Ld^ — dL^, we obtain the following: 

Theorem 2 (Stokes' theorem for nonsmooth chains). If A E X n ,\ and u> E (j^-n+i 
is an (n — l)-form defined on R m , then J A dui = f dA u. 

Historical comments. In 1982 the author found a C 2+a counterexample to the 
Seifert Conjecture [2, 3]. At the heart of her construction is a diffcomorphism of the 
two-sphere with a fractal equator 7 as an invariant set. The Hausdorff dimension 
of 7 is precisely 1 + a. Coarser relations between the differentiability class of a 
function and the topological dimension of a related set were well known in the 
theories of Sard [13], Denjoy [1], and Whitney [15]. Subsequently, these too were 
found to have fractal versions. (See [6, 11].) In each of these theories appears 
the same phenomenon — the smoother the function, the rougher an associated set. 
This paper is a result of the author's search for a general principle underlying this 
duality between dimension and differentiability. 

Fractals are rife in many fields: geometric measure theory, dynamical systems, 
PDEs, function theory, foliations. . . . Except for various dimension theories, how- 
ever, there are few available techniques for investigating them. Usually there is 
a large measure of geometric control over their formation that can be proved or 
reasonably postulated. In at least some interesting cases this ought to allow inte- 
gration of sufficiently smooth forms over the fractals. For example, it would be very 
useful, and perhaps not too difficult, to show that certain solenoid-like invariant 
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sets for flows are domains for integration of 1-forms and that certain closed unions 
of leaves of a foliation by fc-manifolds are domains for integration of fc-forms. Such 
sets could then be proved fairly directly to carry homology cases. Currently such 
proofs (due in various contexts to S. Schwartzmann, J. Plante, D. Sullivan, and 
others) are indirect and hard to come by. It may well be that these new discover- 
ies in geometric integration theory will provide simple but powerful new tools to 
sharpen differentiability results and reduce technicalities in known proofs. 

2. Norms on chains 

We take Whitney's approach in [14] with simplicial n-chains. It provides an al- 
gebraic way of equating chains with common simplicial subdivisions. An n-simplex 
a = po ■ ■ -p n in R" is the convex hull of n + 1 vertices pi in R™, i = 0, . . . ,n. If 
the n vectors {p — Pi} are linearly dependent, then a is degenerate. The order 
of the vertices pi of a nondcgencratc simplex a determines an orientation on a. 
The simplex —a is identified with the simplex with the same pointset as a but the 
opposite orientation. Define a simplicial n-chain A in R n as an equivalence class of 
formal sums aj<7j with real coefficients and oriented simplices in R™ as follows: A 
formal sum ^ aiUi defines a function A: R" — > R by A{x) — ^ ±ai for x e int(ai) 
where +ai is used if <7j is oriented in the same way as R" and — ai is used otherwise. 
Set A{x) — if x docs not lie in the interior of any simplex cr,. Two formal sums 
of oriented simplices arc equivalent if the functions defined by them are identical 
except in sets with n-dimensional Lebesgue measure zero. (Degenerate simplices 
can be used but are equated with the zero n-chain.) The definition of a simplicial 
n-chain A in R™ is independent of the orientation chosen for R™. 

More generally we can define a simplicial n-chain in R m when m> n. For each 
affine n-plane P in R m let C n (P) denote the linear space of simplicial n-chains in 
P. A simplicial n-chain A in R m is an element of the direct sum of the C n (P). 
Let A and B be two simplicial n-chains in R m . If the corresponding summands of 
A and B are equivalent, we write A ~ B. 

Let M n {(j) denote the n-dimensional mass or volume of an n-simplex a. Let 
A be a simplicial n-chain in R™ and n the orthogonal projection onto an affine 
subspace of R m . For < A < n define the projected X-mass of A to be 

M\,it{A) = inf {Y, \ai\(M n (7rai))^ n : ^a 4 a, ~ a} . 

If 7r is the identity, we write M\(A) — M\^{A) for simplicity of notation. Then 
M n (A) is the n-mass of A. It is a norm on simplicial n-chains in R m . 

Integration of smooth integrands over smooth domains. The mass norm 
M n is often used to estimate the integral of a continuous form over a domain with 
finite mass. For example, let A be a simplicial n-chain in R m and lo a continuous 
n-form. Clearly, 



<M n (A)\\u\\ . 



Whitney used a smaller norm on chains to find better estimates for simplicial n- 
chains with large mass. 

Whitney's flat norm. Let A be a simplicial n-chain in R m . Whitney defined the 
flat norm of A as 

\A\ b = mi{M n (A - dC) + M n+1 (C)}, 



4 



JENNY HARRISON 



where the infimum is taken over all simplicial (n + l)-chains C. For example, if A is 
a Jordan curve in R 2 , is bounded by the area of the Jordan domain spanning 
A. If A is an arc in R 2 , consider any simplicial 1-chain B so that A + B is a Jordan 
curve. is bounded by the sum of the arc length of B plus the area of the Jordan 
domain spanning A + B. 

The fiat norm is used to complete the space of simplicial n-chains in R m . De- 
note the resulting Banach space of flat n-chains by X n and the dual space of 
flat n-cochains by X n . Whitney showed that the boundary operator d defined 
on simplicial n-chains is continuous with respect to the flat norm and extends to 
d : X n — > X„_i. 

Let A be a simplicial n-chain and C a simplicial (n + l)-chain in R m . Let u be 
an n-form of class C 1 defined on R m . Whitney used Stokes' theorem for simplicial 
n-chains to show 



Therefore, IJ^^I < |A| b ||w||i. This estimate enabled Whitney to integrate over 
elements of X n . Whitney believed he could generalize Stokes' theorem and wrote 
at the end of his first paper of geometric integration theory [16], "With the help of 
methods described above a very general form of Stokes' theorem may be proved. 
We shall not give details here." However, it turned out that his integral does not 
extend to the Lebesgue integral. Theorem 1 generalizes Whitney's theorem and 
leads to strict generalizations of the Lebesgue integral and Stokes' theorem. 

The natural norms [5]. Let 7 be a Jordan curve in R 3 that has Hausdorff di- 
mension greater than two. This curve not only has infinite arc length but has no 
spanning surface with finite area. Thus a sufficiently close simplicial 1-chain A will 
only bound 2-chains with huge area. Neither the arc length norm Mi(A) nor the 
flat norm \A\^ will give good estimates for integrals over A. The natural norms 
defined below give sharper estimates. 

Since dd = 0, only simplicial n-chains without boundary have spanning (n + 1)- 
chains. Yet all simplicial (n+ l)-chains C can be viewed as "partial" spanning sets 
for an arbitrary simplicial n-chain A. The natural norms take into account the 
(n + l)-mass of these partial spanning sets C and the projected n-mass of what is 
left over, namely, A — dC. (See Figure 1 on page 240.) 

List the n-dimensional coordinate planes of R m from 1 to N — (™) . Let 7Tj be 
the projection onto the ith coordinate plane in this list. 

Let A be an n-chain in R m . If < A < n, define the X-natural norm of A by 



Then \A\\ is proportional to the n-mass M n (A). 

For n < A < m, A e R, define the A-natural norm of A by recursion on the 
integer part of A — n: 




iV 



\A\\ = Y,M X ^{A). 



N 



\A\{ = ^^{Mn,,M-dC) + \C\{}. 



i=l 
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Figure 1 



Each term is infimized over all simplicial (n + l)-chains C in R m . The flat norm 
corresponds to A = n + 1. It is interesting to observe that A may be an integer 
much larger than n + 1 . This allows us to work with domains that are curves in 
M 3 , for example. The A-natural norm depends on the parameter A, the topological 
dimension n of A, and the minimal ambient space dimension m of A. However, if 
A is an n-chain in MP and in M. p , the norms defined are identical. For simplicity of 
notation the dependency on n and m are suppressed. The definition is independent 
of choice of coordinates. That is, under a C x ~ n+1 change of coordinates with 
bounded derivative, the A-natural norm changes to an equivalent norm. 

As does Whitney, we complete and obtain a Banach space X n ,\ of X-natural 
n-chains and its dual space X n,x of X-natural n-cochains. 

Remark. There is a one-parameter family of homology and cohomology groups 
associated with the A-natural norm, H n ,\ and H n ' X . 

Now let w be a 1-form defined on R m of class C 1 and A G Xi t \. We have seen 
that the integral of u over A can be estimated using either arc length or the flat 
norm, both of which may be very large. Theorem 1 shows that if u is of class C 2 , 
the integral can also be estimated using the 3-natural norm which may be much 
smaller. 



To prove Theorem 1, we iterate a modified version of the Stokes' argument as 
seen in (*) (see [5] for details; also see the example below). The reader should take 
note that Stokes' theorem, as such, may be applied only once since ddco = and 
ddA = 0. However, the components of dui may not be exact, and we may apply 
the exterior derivative to them. This means we may continue to use the exterior 
derivative and iterate the "Stokes' process" in a nontrivial fashion for as many times 
as lu is diffcrcntiable. We also use partial spanning sets of spanning sets. 

Typically, an integral J A to is treated with two basic methods: If A = dB and 
Lu is smooth, one may "go forward" and integrate dio over B. Alternatively, if the 
form to = du, one may "go backward" and integrate v over dA. Sometimes neither 
of the equivalent new integrals is easier to calculate. Our methods allow one to 
go forward many times under suitable conditions, taking the exterior derivatives of 
forms in the hopes of finding an equivalent integral that may be calculated. It is 
worth mentioning that it is also possible, under suitable conditions, to go backward, 
taking antiderivatives of forms many times, and find an equivalent integral that may 
be calculated. 



(**) 
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Example. There exists a self-similar Jordan curve 7 in the unit cube with Haus- 
dorff dimension > 2 so that each projection 7^(7), i — 1,2,3, onto the ith two- 
dimensional coordinate plane is an immersed curve bounding an infinite 2-chain 
with finite 2-mass < 1. (The curve 7 is constructed by adding homothetic replicas 
of the curve a to itself in Figure 1 so that 7 is embedded and has Hausdorff di- 
mension ~ log 11/ log 3 > 2.) A simplicial 1-chain A can be found approximating 7 
along with simplicial 2-chains d spanning A, i = 1,2,3, such that M 2j7ri (Ci) < 1. 
(Each surface Ci is a sum of homothetic replicas of the surfaces a and r in Figure 
1, chosen to minimize the projected 2-mass of C,.) Furthermore, the simplicial 
3-chains Dij spanning C\ + Cj have M^(Dij) < 1. Even though the norms Mi (A) 
and \A\ b are both large, the 3-natural norm \A\\ is not, for 

3 3 
|^| 3 <^(M 2;7r .(C j ) + |A j | 3 )<12 and thus |^<^|C^<36. 
j=i »=i 

This curve 7 cannot be a domain for the classical Lebesgue integral and cannot be 
treated using Whitney's methods. We show in §3 how to choose Ak approximating 
7 and define J lo = lim J A ^ u). 

We verify that Theorem 1 (sec (**)) is valid for this example. It suffices to 
estimate | J A LVi\ = | J c duJi\ where u>i is a component of ui. Note that Ci + Cj — 
dDa. Thus 



/ duji 


3 


/ 


+ 


/ d(dLUi)j 


<\A\l\H\ 2 


JCi 




Jc 3 




J 





Relation to the Lebesgue integral. Our integral of smooth forms contains as a 
special case the Lebesgue integral: First associate an L 1 function / : M 1 — > M 1 with 
an element of Xi_2 as follows. Assume for simplicity that / is the characteristic 
function of a bounded, measurable set E. Let S denote the region under the graph 
T of / and above the x-axis, oriented positively. Let m„ denote n-dimensional 
Lebesgue measure. Since / G L 1 , m 2 {S) < 00. For each k there exists a union Pk 
of disjoint intervals such that E C Pk and mi(Pfc) < 00. Furthermore, f] Pk = E 
except for a zero set. Subdivide each of the rectangles of Pk x / into finitely many 
2-simplices and orient each positively. Denote the formal sum of these oriented 
2-simplices by Qk- Since m,2{Pk x 7) < 00, each Qk is an element of X 2 ,2- Since 
rri2(S) < 00, the sequence {Qk} is Cauchy and thus converges to Af e A 2 ,2- It can 
be shown that dAf e X\ t 2 is canonically associated with /. Now it is well known 
(see Saks [12]) that the Lebesgue integrals J R / and J s dxdy are equal. Since S is 
canonically represented by Af G X2.2, the Lebesgue integral J s dxdy equals our 
integral J A ^ dxdy. Since ydx is analytic and dAf e Xi^, Theorem 2 implies the 
Lebesgue integral J R / equals our integral J r y dx. (This example can be treated 
with the flat norm alone.) 

3. Domains of integration 

Which oriented n-dimensional topological submanifolds M C K m with boundary 
can serve as domains for integrals of smooth forms w of class C A ~™, n < A < m? 
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The idea is to choose appropriate simplicial approximators Ak for M that converge 
in X n ^\ for A sufficiently large and define f M w — lim f Ak ^. (See [4] for some 
numerical methods.) But choosing Ak takes care. For example, let Mi be a Jordan 
curve in M 2 with positive Lebesgue area. Any Cauchy sequence of simplicial 1- 
chains inside Mi will have different limit point in from a Cauchy sequence of 
1-chains outside M\. The area of Mi itself contributes an unavoidable error. 

Consider a compact oriented n-manifold M C R m with boundary. There is 
an algorithm for constructing a particular sequence of simplicial n-chains in the 
nerves of coverings of M by boxes in cubic lattices, the inverse limit of these chains 
representing the fundamental Cech homology class of (M, dM) . If this sequence of 
binary approximators converges to some ?/>(M) e X n ^\(R m ), then M is a (A,n)- 
set. The algorithm commutes with the boundary operator on chains. Examples of 
(A, n)-sets include some planar Jordan curves with positive Lebesgue area and the 
above example of a Jordan curve in 3-space with Hausdorff dimension > 2. If M is 
a (A, n)-set and u is an n-form of class C A ~™, we can define J M u = f^ M ^ oj and 
apply Theorem 2 to conclude: 

Theorem 3 (Stokes' theorem for (A, n)-sets). If M is a (A, n)-set andu; e C x ~ n+1 , 



Theorem 3 extends the classic Green's theorem for codimension one boundaries 
7 and C 1 forms u>. Holder versions of the classic Green's theorem relating the box 
dimension of 7 to the Holder exponent of to are proved in [7, 8, 10]. 

There are many examples of (A, n)-sets: 

(i) There exist Jordan curves in R n with Hausdorff dimension d that are (d, 1)- 
sets for every 1 < d < n. 

(ii) Every hypersurface immersed in R m with Hausdorff dimension < to is an 
(to, to — l)-set. 

(iii) Compact oriented n-manifolds immersed in R m that are locally Lipschitz 
graphs are (n, n)-sets. 

(iv) Compact oriented 1-manifolds immersed in M 2 that are locally graphs of 
Holder functions with exponent a arc (2/(1 + a), l)-sets. 

Remark. Compactness in examples (ii), (iii), and (iv) can be relaxed. An arc in M 2 
that spirals to the origin passing through the x-axis at x n = 1/y/n, n > 1, is not 
a (A, l)-set for any 1 < A < 2. Other domains of integration include boundaries of 
open sets such as the topologist's since circle and the Denjoy Cantor set. 

Integration on manifolds. We pass to smooth manifolds by way of local coordi- 
nates, since the A-natural norm changes to an equivalent norm under a (A — n + 1)- 
smooth change of coordinates with bounded derivative. One must also extend the 
definition of these norms to chains in proper open subsets of R m as Whitney did 
for his flat norm. The integral can then be defined for smooth forms defined on 
smooth compact m-manifolds M over domains that are "rough" subsets of M. 
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